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First principles density functional calculations of the electronic and magnetic properties of spinel- 
structure LiV2 04 have been performed using the full potential linearized augmented planewave 
method. The calculations show that the electronic structure near the Fermi energy consists of a 
manifold of 12 bands derived from V t2 g states, weakly hybridized with O p states. While the total 
width of this active manifold is approximately 2 eV, it may be roughly decomposed into two groups: 
high velocity bands and flatter bands, although these mix in density functional calculations. The 
flat bands, which are the more atomic-like lead to a high density of states and magnetic instability of 
local moment character. The value of the on-site exchange energy is sensitive to the exact exchange 
correlation parameterization used in the calculations, but is much larger than the interaction between 
neighboring spins, reflecting the weak coupling of the magnetic system with the high velocity bands. 
A scenario for the observed heavy fermion behavior is discussed in which conduction electrons in the 
dispersive bands are weakly scattered by local moments associated with strongly correlated electrons 
in the heavy bands. This is analogous to that in conventional Kondo type heavy fermions, but is 
unusual in that both the local moments and conduction electrons come from the same d-manifold. 



I. INTRODUCTION 

Heavy fermion (HF) materials are typically inter- 
metallic compounds containing Ce, U or Yb atoms. They 
are characterized by the usual Landau Fermi liquid scal- 
ing properties, but only at very low temperature (be- 
ginning as low as 0.3 K, depending on the material) 
and with extraordinarily strongly renormalized effective 
masses, m* w 100— 1000m e .Ertj They show apparent local 
moment paramagnetic behavior with strongly increasing 
spin susceptibility, x, and specific heat coefficient 7 with 
decreasing temperature, but do not order magnetically 
and eventually settle at low T into a state with constant 
7 and x an d Wilson ratio near unity. Following exten- 
sive investigation over many years, a basic understand- 
ing of the phenomena has been established. The origin 
is a many body effect associated with the interaction of 
itinerant conduction electrons with strongly correlated f- 
electrons on local moment rare earth ions. Thus the dis- 
covery by Kondo and co-workers of HF behavioiu with 
7 » 0.42 J/mol K 2 , in the transition metal oxide L1V2O4 
was both remarkable and unexpected as the material does 
not appear to fit into this framework. 

As may be expected, this discovery has led to a series 
of detailed experiments that have confirmed the origi- 
nal result with properties-characteristic of the heaviest 
f-electron HF materials JjTLTJ although it should be noted 
that Fujiwara and co-workers proposed an alternate ex- 
planation of their NMR data within a spin-fluctuation 
framework.Q 

Two first principles band structure calculations have 



recently appearedLM, though both use spherical approx- 
imations for the potential. These reasonably agree with 
each other but not as well with the full potential calcula- 
tions presented here. Their conclusions are, however, dif- 
ferent: Eyert et alJld suggest that L1V2O4 is not a "true" 
HF material in the sense that they ascribe the specific 
heat enhancement to spin fluctuations arising from frus- 
trated-antiferromagnetism on the V sublattice; Anisimov 
et al.t 2 ] advocate a separation of V d electrons into two 
sets, a localized and an itinerant one, playing the role of 
the f and conduction electrons in convention&LHF mate- 
rials. A related model is discussed by Varma.Ej We argue 
that the true HF picture is likely, since it is compatible 
with the electronic structure obtained and is consistent 
with the experimental situation. 

L1V2O4 occurs in an undistorted cubic fee spinel struc- 
ture as shown in Fig. 0. The V+ 3 5 (3d 1.5) 10ns are 
located in very slightly slightly distorted O octahedra. 
These V ions form corner sharing tetrahedra ordered on 
the fee lattice such that each V has six neighbors of the 
same kind. Thus the V sublattice can be viewed as 
an fee lattice with every second atom removed. Alter- 
nately, it can be viewed as a collection of chains running 
along all the (110) type directions with each V partic- 
ipating in three chains. This structure leads to strong 
geometric frustration of antiferromagnetic (AF) interac- 
tions, although of course ferromagnetic (FM) interac- 
tions are unfrustrated. As noted, L1V2O4 is paramag- 
netic at— low temperature. However, Krimmel and co- 
workersEj have reported quasi-elastic neutron scattering 
measurements that show predominantly FM fluctuations 
above 40K with substantial observable AFM fluctuations 
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only at lower temperatures, where they significantly con- 
tribute to the spectra. 



II. BAND STRUCTURE AND MAGNETISM 

In this structural configuration, we find that the O 
2p bands are below and well separated from the transi- 
tion metal derived bands so that both the Li-0 and V-0 
bonds are strongly ionic. The higher lying 3d derived 
orbitals are crystal field split into a 6-fold degenerate (3 
per spin) ti g manifold and a higher lying 4-fold degen- 
erate e g manifold. Neither the crystal structure nor the 
electronic configuration is particularly favorable for sub- 
stantial V 3d-0 2p hybridization, so narrow d bands are 
obtained with clean crystal field gaps between the ti g and 
e g manifolds. We do not get noticeable mixing of these 
manifolds. Thus the electronic structure relevant to low 
energy excitations may be described as a manifold of 12 
narrow ti g bands 1/4 filled with 6 electrons (note that 
there are 4 V ions per cell) with no other nearby bands. 
Much of the bandwidth is due to direct V-V hopping. 
At first glance it is hard to recognize the connection be- 
tween this electronic structure and that of the f-electron 
HF compounds. 




FIG. 1. Structure of LiV204 showing the local coor- 
dination of a V atom. Note that central V in the figure 
is at the center of an O octahedron oriented along the 
100 axes and on the joint vertex of two corner sharing V 
tetrahedra, leading to the rhombohedral site symmetry. 

We find, using density functional calculations with the 
full potential linearized augmented planewave (LAPW) 
method,t-irtZI that the ti g manifold in fact contains two 
different types of carriers as was also emphasized by 



Anisimov et al.EH. Calculations with the tight-Jainding 
linear muffin tin orbitals (TB-LMTO) methodrJ were 
used to analyze the band symmetries. One carrier type 
is analogous to the itinerant conduction electrons of con- 
ventional HF materials and the others can play the role 
of the strongly correlated local moment f-electrons. This 
provides a way of obtaining conventional HF behavior 
in LiV2 04 albeit in an unconventional way, at least as 
regards the nature of the bare states. 




FIG. 2. Calculated band structure 

of non-spin-polarized LiV2 04. The Fermi energy is de- 
noted by the dashed horizontal line at 0. 

The calculated local spin density approximation 
(LSD A) band structure for non-spin-polarized (NSP) 
LiV204 is shown in Fig. ^. The corresponding elec- 
tronic density of states is given in Fig. [| The O 2p 
bands span the range from approximately —8 to —3 eV 
(relative to Ep ). The 12 i 2 g bands are separated from 
these by a gap and have a width of slightly less than 
2.5 eV. The e g manifold is quite narrow, characteristic 
of the bonding topology of spinels, which features bent 
V-O-V bonds, unfavorable for band formation via e g —pa 
hybridization. These are again well separated from the 
active ti g bands by a gap of slightly less than 1.5 eV. The 
dispersion of the ti g bands is in fact mostly derived from 
V-V hopping as was verified by removing the O orbitals 
from TB-LMTO calculations. A blow-up of the band 
structure showing the t2 g bands, which dominate the low 
energy physics, is given in Fig. ^. Although a high on- 
site Coulomb repulsion, U , is not expected in early 3d 
compounds like this, significant correlation effects on this 
band structure should not be excluded since the .active 
manifold is narrow. Photoemission measurementsE2l im- 
ply U as 2 eV, which is comparable to the bandwidth, 
and therefore consistent with at least moderate correla- 
tion effects in this multi-band system. 

Spin polarized calculations, both within the LSBA and 
with generalized gradient approximations (GGA)cJ show 
a considerable magnetic instability, of strong local mo- 
ment character. Calculations were done for ferromag- 
netic, antiferromagnetic (AF — 2 of 4 spins flipped per 
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unit cell) and ferri-magnetic (FiM — 1 of 4 spins nipped) 
configurations. Very low Hellman-Feynman forces on the 
symmetry unconstrained O coordinates were obtained in 
the experimental crystal structure, used in these calcula- 
tions, independent of the magnetic order, supporting the 
experimental structure and implying only weak magneto- 
elastic effects as may have been anticipated based on 
the weakly hybridized electronic structure. Within the 
LSDA, a FM spin polarization of 1.15 /^b/V with an en- 
ergy 60 meV/V below the NSP energy is obtained. The 
AF and FiM energies were both 5 meV/V lower than the 
FM. Thisas qualitatively different from the results of Ey- 
ert et alSJi, who found energies and moments closer to 
itinerant magnetic system consistent with strong scatter- 
ing of carriers by spin fluctuations at the LSDA level. 




FIG. 3. Electronic density of states of non-spin -po- 
larized LiV2C>4 in states per eV per unit cell (4 V atoms) 
as a function of energy in eV. The Fermi energy is at 0. 



of states. This gave a FM magnetization of slightly more 
than 0.5 and slightly more than 1 /Xb/V respectively, in- 
dicating that the changes in the bare NSP bands are the 
more important factor, although it should be noted that 
with the full exchange splitting the bands are non-rigid 
so the extended Stoner is not fully valid and the result 
should be regarded as qualitative. 
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FIG. 4. Calculated band structure of non -spin -po- 
larized LiV204 showing the t2 9 manifold. The sizes of the 
points denote the calculated relative amounts of a% g char- 
acter in the states. The Fermi energy is denoted by the 
dashed horizontal line at 0. 



The GGA results are qualitatively similar to the 
LSDA, but are somewhat more magnetic. In the GGA 
the FM spin polarization is increased to 1.4 /i_b/V, which 
is a greater than usual sensitivity to the exchange cor- 
relation functional. The magnetic properties are also 
unusually sensitive to other details of the calculations: 
for instance, the atomic sphere approximation (ASA), 
used in Refs. [y],[l^, as well as in our LMTO calculations, 
leads to substantial underestimate of the magnetic mo- 
ment, namely 0.5 /xb/V. This is also unusual and calls 
for caution when using ASA calculations for quantita- 
tive estimates of the electronic parameters (cf. Ref. 12). 
This is related to the physics discussed below. First of 
all, the magnetization comes almost entirely from one 
d sub-manifold, resulting in a highly non-sphcrical ex- 
change potential, which is spherically averaged in non- 
full-potential calculations reducing its effect. Secondly, 
the NSP energy bands differ such that the full-potential 
band structure has higher values of the density of states 
around the Fermi level, and therefore a stronger Stoner 
instability. To test this, extended Stoner calculations 
were performed with the same Stoner / (from the TB- 
LMTO), but the LSDA TB-LMTO and LAPW densities 



The calculated spin density (Fig. g) is remarkable, in 
that despite the nearly cubic environment, it strongly re- 
flects the weak rhombohedral symmetry. This is a result 
of the fact that the hopping that gives rise to the t^g 
dispersions is largely V-V and that the V neighbors are 
a rhombohedral environment. The primary pseudocubic 
quantization axis is along the Cartesian 100 directions 
because of the octahedral mostly O derived crystal field 
that splits the ti g and e g manifolds. In the rhombohedral 
site symmetry there is a further quantization of the ti g 
space based on the rhombohedral 111 axis, into a double 
degenerate e g space and a single degenerate a\ g mani- 
fold. In order to distinguish this e g from the higher lying 
e g manifold split off by the octahedral crystal field, we 
will denote it e' g in the following. The e' g and a\ g are 
at the same energy in LiV204 but are locally orthogonal 
spaces. The remarkable result, evident in Fig. [| is that 
the magnetization is derived almost entirely from the a\ g 
sub-space. 
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FIG. 5. Calculated spin density of LiV204 with a fer- 
romagnetic ordering. The top panel shows an 0.1 fiB /a. it. 3 
isosurface viewed along a 111 direction, the middle is the 
same view blown up around a single V ion, while the bot- 
tom panel is the same on a -110. Note the predominant 
but not pure ai g character. The small e' g admixture is 
responsible for the triangular rather than full symmetry 
viewed from the top (middle panel) and the buckling of 
the central "doughnut" in the side view (bottom panel). 



We obtain almost complete polarization of the relevant 
ax g sub-manifold in our full-potential LSDA and GGA 
calculations. Besides the higher resulting magnetic mo- 
ments, we obtain much clearer strong local moment be- 
havior in these calculations. The V moments as charac- 
terized by integrated spin magnetization in the V LAPW 
spheres for the AF and FiM cases are within 9 percent of 
the FM value. The FiM and AF cases were energetically 
degenerate to within the accuracy of the calculations at 
5 meV/V below the FM energy, consistent with weak 
AF exchange interactions and magnetism that, as noted 
above, is very strongly local moment in character. 

As mentioned, the dispersion of the t2 g bands is nearly 
exclusively derived from the V-V interaction. This in- 
cludes the nearest neighbor dda and ddn hoppings. In- 
verting the TB-LMTO Hamiltonian we get a rough esti- 
mate, t dd(7 /t dd7T ~ 3.5 (and t ddS < t ddrr ). The 12x12 TB 
Hamiltonian in this case is defined by the overlap inte- 
grals for each V-V bond, e.g., for the (110) bond t xy - xy — 
jt dda , t yz - xz = t ddlTl t xy - xz = 0. An interesting prop- 
erty of this Hamiltonian is that when t dda — %t ddlx , two 
bands are flat. These bands are nonbonding combina- 
tions of e' g states of the four V ions. In actual calcula- 
tions this condition is satisfied only very approximately 
(within k 40 percent), but the corresponding states are 
still very flat (the width is 0.2 eV) and located about 0.8 
eV above the Fermi level. Full-potential LAPW calcula- 
tions reveal also two more relatively flat bands of the e' g 
character right below these, but still about 0.7 eV above 
the Fermi level. On the other hand, a closer look reveals 
four more heavy bands (width ss 0.5 eV) that are located 
around the Fermi level. Those have predominantly a\ g 
character and give rise to the high concentration of a± g 
character around Ep and the high density of states in the 
same region with the resulting magnetic instability (cf. 
Fig. ||). The remaining bands are mixed in the LSDA; 
e' g character dominates and they have greater dispersion. 



III. SCENARIO FOR HEAVY FERMION 
BEHAVIOR 

The variation of the effective mass is substantial: near 
the bottom of the ti g manifold it is close to the free elec- 
tron mass (which is unusually small for d-bands), while 
near the Fermi level it is 4-5. From the transport point 
of view, one can speak about two groups of carriers: four 
light bands (comparable with the sp bands in transition 
metals), and four heavier bands (with their width of « 0.5 
eV). Additionally, there are four flat bands at about 0.7- 
0.8 eV above Ep. The bands providing the heavy carriers 
like the light bands (width 2 eV) cross the Fermi level but 
are responsible for the calculated local moment magnetic 
instability. These heavy bands near Ep are dominated 
by aig character, while the remainder of the bands, al- 
though mixed are dominated by e' g character. 

In the scenario proposed here, the states relevant to the 
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transport properties are from the heavy and light bands 
near Ep] these may constitute the two classes of carriers 
mentioned above. This connects with conventional HF 
materials, in that LSDA calculations for these materials 
also show band structures consisting of broader conduc- 
tion bands and narrow (W <C U) bands at Ep associated 
with local moments. In conventional HF materials these 
two classes of bands are coupled, but weakly. The first 
evidence for this in L1V2O4 is through the small (but 
no doubt overestimated in LSDA) inter-atomic exchange 
interactions as seen by the near degeneracy of the var- 
ious magnetic configurations studied and the fact that 
although the ferromagnetic band structure has exchange 
splittings in all bands, very little polarization of the e' g or- 
bitals is induced as discussed above. Further, more clear 
signatures emerge from the transport functions, particu- 
larly N(Ep)(vp) oc Vujp oc V(n/m) where Vp is the band 
velocity on the Fermi surface, V is the unit cell volume, 
Lu p is the plasma frequency, n/m is the optical n/m and 
the average () is over the Fermi surface. This function is 
dominated by more dispersive bands due to the (vp) fac- 
tor, so that changes with magnetic order are reflective of 
the interaction between the low dispersion states respon- 
sible for the magnetism and the more dispersive conduc- 
tion aig derived bands. We obtain N(Ep) of 7.26 eV" 1 
on a per formula unit basis (2 V atoms) for the param- 
agnetic calculation, and 7.23 eV^ 1 for the FM majority 
spin channel, 1.58 eV _1 for the FM minority channel, 
2.96 eV _1 for the FiM majority spin channel and 7.24 
cV _1 for the FiM minority channel. The correspond- 
ing Fermi velocities (cm/s) are 0.77 x 10 7 , 0.30 x 10 7 , 
1.37 x 10 7 , 0.60 x 10 7 and 0.75 x 10 7 . The LSDA values 
of N(Ep){v F ) stay within 20 percent of the paramagnetic 
value. 



So far all our results emerge from first principles den- 
sity functional calculations primarily within the LSDA 
using the full potential LAPW method. This includes 
the separation into two different types of carriers, a\ g like 
and e' g like and the strongly spin polarized a\ g like local 
moment magnetism. Coulomb correlations as included 
approximately in the LSDA+U method are not required. 
This is of some significance because in the LMTO calcu- 
lations of Ref. [l^ a U of 3 eV was used to obtain these, 
and this U is about 50 percent greater than the value 
implied by photoemission and which is also a value that 
would split the e' g space into strongly correlated upper 
and lower Hubbard bands. On the other hand, with the 
screened U less than or equal to the value indicated by 
photoemission, as is normal, i.e. U < 2 eV, the light e' g 
derived bands would be weakly to moderately correlated 
and metallic as in conventional HF while only the heavy 
aig bands would be strongly correlated. In fact, the ac- 
tual situation may be more complicated because the two 
relevant types of carriers are derived from the same d 
manifold and therefore may be more strongly coupled 
than in conventional HF, as discussed below. 



IV. DISCUSSION 

Even in weakly correlated materials, LSDA calcula- 
tions overestimate the extent of hybridization between 
different orbitals. Thus one may anticipate the follow- 
ing artifacts in the LSDA calculations: (1) the heavy 
aig bands could be narrower than calculated and (2) the 
mixture between a\ g and e' g bands could be weaker than 
calculated and (3) as a result the dominance of a\ g char- 
acter in the spin density could be stronger than that ob- 
tained from the LSDA. This is like the corresponding ar- 
tifacts in band calculations for f-band metals where the 
hybridization of the f-bands with the conduction bands 
is overestimated in the LSDA. These are partially lo- 
cal effects related to LSDA errors in the description of 
the bonding and not directly to Hubbard type correla- 
tions controlled by U/W (the same effects are present 
in wide band materials, where, for example, LSDA es- 
timates of interatomic magnetic couplings, J are often 
overestimated) . Additionally, the effect of Hubbard type 
correlations, which are neglected in band structure calcu- 
lations, may be expected to be strong for the heavy mass 
carriers as taking LL< 2 eV, as indicated by photoemis- 
sion measurementst£l, yields U/W < 4. At such values of 
U/W one may expect a substantial suppression of charge 
fluctuations involving the heavy mass bands along the 
lines of the normal f-electron HF systems. However, the 
real situation is undoubtably considerably more complex, 
as some Hubbard correlation effects may be present in the 
bare light carrier bands where U/W < 1 and exception- 
ally strong off-diagonal Hubbard effects may be expected 
as well since the relevant bands (heavy and light) are 
both derived from the V t2 g space. It will be quite inter- 
esting to construct and study many body models of this 
type in the HF regime. Only the t2 g orbitals and the in- 
teractions among them should be needed. There should 
be four light bands and four related flat bands away from 
Ep per unit cell. These should be formed from two or- 
bitals per V using V-V hopping as described in the tight 
binding discussion above. The light bands, which provide 
the carriers, should have width 2 eV, U/W < 1 and con- 
taining w 0.5 electrons per V. In addition, there should 
be a localized set of atomic-like orbitals (one per V) with 
occupation le and bare band width « 0.5 eV. Addition- 
ally, the off-diagonal Coulomb repulsion coupling these 
two groups of states should be comparable, i.e. » 2 eV. 
LiV204 can be made quite cleanly and is well suited to 
detailed experimental investigation, and this may make 
it an especially interesting system for experimental and 
theoretical study. 
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